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Abstract 

^ ' We relate possible cosmological variations in the mass ratio /i = mp/me 

, and the fine structure constant a to long-range composition-dependent 

^^O ' forces mediated by a scalar field. The differential acceleration 77 in Eotvos- 

- type experiments is bounded below by 10~^^, except in cases where one or 

more scalar couplings vanish. We consider what values for these couplings 
, could arise from unified theories. By considering the contribution of the 

' scalar field to the cosmological energy density we use bounds on r] to put 

upper bounds on the current rate of change of /l( and a. 

p. 



O-i. 1 Introduction 



The constancy of the coupling strengths and particle masses in the Standard 
Model Lagrangian is an assumption that should be tested. If a variation existed 
^ ' it would violate the principle of Local Position Invariance contained in the 

■ Einstein equivalence principle. Nevertheless consistent, relativistically covariant 

theories with varying couplings can be written down by including scalar degrees 
of freedom, in which the variation arises due to the cosmological solution for 
the scalar (s). See for example US HZl HH EH] • 

There are also many situations in cosmology where gravitational effects 
appear to behave differenty from what a naive application of GR would produce, 
which may also motivate the existence of cosmological scalar fields. The first 
is inflation (see for example ^01)) which seems to require a scalar field with a 
peculiar kind of potential. Closely related, though happening on an immensely 
slower timescale, is the current apparent accelerated expansion of the Universe 
^DEIl) which may also be related to scalar field dynamics [T^ 1^. 

The violation of the Einstein equivalence principle makes itself via scalar 
mediated forces, thus gravity becomes a generalised scalar-tensor theory. So far, 
precision tests of gravity have revealed no deviation from General Relativity |311 
I22j . in particular concerning the question we will discuss here, the equivalence 
of gravitational and inertial mass (universality of free fall) 'T^ . The motivation 
for continuing to consider deviations from GR then comes into question. 

* email: dent@thphys.uni-heidelberg.de 
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The current motivation comes from astrophysical measurements of various 
coupling strengths and mass ratios deduced from atomic and molecular absorp- 
tion (or emission) wavelengths. These probe back to redshifts up to about 4 and 
have now reached the sensitivity of 10~^ or better; moreover some measurements 
of the fine structure constant a 1 , and most recently the proton-electron mass 
ratio n = mp/me show potentially significant deviations. The fine structure 
constant has been the subject of debate, with other recent studies obtaining 
negative results [SI E] , some at a precision greater than or equal to previous 
evidence for variation. 

Bounds on variation of a from nuclear physics effects [HI [Hj have also been 
derived; however they are subject to considerable uncertainty if one also allows 
other quantities such as quark masses to vary (see for example j29j). 

It was pointed out in J3] that a fractional variation of a on the order of 
10~^ required a scalar coupling to electromagnetic energy large enough such 
that the differential acceleration of two bodies in the gravitational field of the 
Sun (including scalar forces) is of order 10~^^, which may be detectable with 
the STEP experiment. A more detailed analysis of scalar field models in jl7j . 
including also the electrostatic energy of nuclei, found a rather more restrictive 
bound from non-universal free fall, generally referred to as Weak Equivalence 
Principle (WEP) violation, although. Other works have related a variation in 
a to scalar field models which describe the current acceleration of the Universe 

[iHi[ini[2ni. 

As pointed out in ^21,, one reason for believing the result of [T^ to be an 
underestimate of WEP violation is that the authors did not include possible 
scalar couplings to other Standard Model sectors which are likely in unified 
models (see also [21 [21I1S1 HZI UHl UHl)- Clearly also a scalar coupling to 
electromagnetism alone cannot produce a variation in /x at the same level as 
that in a. 

Hence in this paper we redo and extend the analysis to include other possible 
types of scalar coupling through which varying a and ^ can be related to WEP 
violation. Current limits on differential gravitational acceleration are at the 
level of 10~^^ - 10~^^ (7j for the parameter rj defined as 

_ |ai - 02! , , 

^ = 2| : 1 (1) 

\ai + 02! 

where two test bodies of different composition have accelerations ai, 02 towards 
a known gravitational source. 

In the first section of the paper we make a preliminary estimate of the size 
of differential acceleration r/ due to (/)-mediated forces, allowing for varying /i, 
and find it may be orders of magnitude greater than 10"^*^, and possibly visible 
with the MICROSCOPE experiment with sensitivity 10~^^. Exceptions occur 
if some scalar couplings vanish exactly. 

In the second section we consider what scalar couplings consistent with vary- 
ing (J, may emerge from unified models both with and without supersymmetry, 
using the detailed analysis of [IHj • We find it unlikely that the scalar coupling 
to any particle will vanish. In fact some unification scenarios imply that WEP 
violation may be at the current limits of detection, or even be ruled out by 
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Eotvos-type experiments. Only a restricted range of unified scenarios can be 
consistent with reported variations of a and n, if these are to arise from a scalar 
field varying over cosmological time. 

In the last section we examine the cosmological evolution of the scalar field 
in more detail and suggest a direct test of unified scenarios by comparing WEP 
violation with atomic clock measurements of the constancy of fundamental pa- 
rameters. 

1.1 Consistency of "varying constant" observations 

Clearly, the experimental situation on possible variation of a is currently unsat- 
isfactory, while that for fi requires further independent testing. We discuss two 
points: first, whether nonzero variations in both a and ^ are inconsistent with 
other measurements of dimensionless fundamental quantities; second, whether 
recent results on a are directly comparable. 

The quantities y = a^Qp and x = a^gpme/mp have been measured with 
good precision using quasar absorption El • Prom two systems at redshifts 
below 0.7 the fractional variation of y is bounded at the level of 5 x 10~^, hence 
is not inconsistent with the claim Aa = (— 0.57ib0.11) x 10~^ (made for systems 
over a wide range in redshift) given constant gp. ^ The result for x indicates 
a slight but not significant increase into the past, Ax/x = (1.2 ib 1.0) x 10~^. 
Again given constant Qp, this is difficult to reconcile with a decrease in a and an 
increase in ^ as claimed in '2^ at the level A/i//x = (2.4 it 0.6) x 10~^. Hence we 
cannot assume that both a and ^ have exactly the claimed nonzero variation. 
Rather we will make estimates assuming that at least one nonzero variation 
exists (usually ^) at the claimed level. 

The most recent negative results for Aa are based on very few absorption 
systems, the improved sensitivity being ascribed to more accurate spectroscopy. 
The Webb group results by contrast involve a sample of about 150 systems.^ 
This makes it possible to compare the estimated error in each systemstatistical 
scatter over the sample, and carry out checks against systematic error. One 
outstanding question is the statistical error due to fitting the velocity profile of 
the absorbing system; see for example the first reference of 1^. It may require 
a large number of systems to be analyzed to settle the question. 

2 Preliminary estimate of j] 

The estimate proceeds along the same lines as ^S], but with some important 
differences. In particular the electron content of matter must be explicitly 
included since the variation of /i involves rUe- 

We posit a light neutral scalar field (j) which is canonically normalised, re- 
defining the origin such that the value of (j) vanishes today, and work in the 
Einstein frame. It will be convenient to use a dimensionless scalar = (p/Mp 

^Possible variation of g-factors was discussed in |30| : any effect on Qp is expected to be 
small. 

^Analysis of more recent observations is ongoing (M. Murphy, private communication). 
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where Mp is the reduced Planck mass (SttG) The scalar couples to matter 
in the low momentum limit via effective operators 

J^<f>m = -ruefl [l + Ae(A] 66 - nipfl [l + Xp(j)] pp - ninfl [l + A„(^] nn (2) 

where m^fi etc., are the present-day values of particle masses, and to electro- 
magnetism via 

C^em = -\[^ + Aem0] F^.F^^" (3) 

thus the fine structure constant is given by a = ao(l + Xem4')- Note that 
Ap and A^ can differ from one another only due to light quark masses and 
electromagnetic effects, i.e. isospin violation. We may reexpress the nucleon 
mass terms as 



ruN [l + Xiff] {pp + nn) + 



rriN 



+ Xv<P 

rriN 



{—pp + nn) (4) 



where m at is an averaged nucleon mass (we drop the "0" suffix except in case 
of ambiguity), 5m„p is the nucleon mass difference and 

_ Xpnip + XnTUn _ -XpUlp + Xnnin 
Ai = , Av = • 

niN nriN 

Note that this form gives only the leading linear dependence on 0. We could be 
more general and replace Ae and other couplings with general functions Xe{(f>) 
(etc.) defined via 

dln{me/Mp) - 
" ^^^'^^ 

allowing a nonlinear dependence, which may be likely if the value of </) changes 
by more than a small fraction. 

For a given test object we have n^ = np electrons and protons and n„ neu- 
trons, with the "proton fraction" fp defined as np/ {np + nn), equal to the mass 
averaged value of Z/A, and "neutron fraction" = 1 — fp. The nuclear elec- 
tromagnetic binding energy also makes a contribution femfnN Z{Z — 1)A^^/^ 
to the mass, where fem — 7 x 10~^. The ratio of the (;/)-mediated acceleration 
of a test body to the gravitational acceleration then 

a0 f fp{Xeme + Xpmp) + f^Xnmn {Z{Z -l)A-^/^)'\ 

J em^e 



ttgr y niN {A}'' J 

( fp{Xeme + XpUlp) + fnXnmn , , (Z(Z - 1)^'^^) ^ 

\ ruN {A) J 

where the superscript s refers to the common source. If, as in some experiments, 
this is the Sun then the electromagnetic contribution to the source coupling is 
negligible. Note that if electrons are in relativistic motion their scalar charge 
is "diluted" by a boost factor, as derived for example in but to obtain 
conservative bounds on scalar couplings we will take the nonrelativistic limit. 
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For two test bodies possessing different fp, Z and A one can determine the 
ratio 



V 



where <C a^r- We then obtain 



a, 



gr 



(6) 



Ae h Ap 

rriN 



+ rx +A f iz{z-i)A^yY .^,^ 

' J n^n ~r ^emj em , I ^ 



rUe nip rrin 
Ae h Ap An 

m^v ■'TiAT 



+ Aem/emAi2 



Ae Ai A^ 

1837 T ~ 



Ae 



1837 



A^, 



+ Aem/emAl2 



{A) J 

(A) 



(7) 



where the notation A12 denotes the difference between the two test bodies and 
we take them to have Z ^ 1. The last hne of Eq. Q follows since we find 
the strongest bound from an experiment using "Earth-like" (mostly iron) and 
"Moon-like" (mostly silica) test bodies accelerating towards the Sun, for which 
~ 1; the test bodies differ by Aufp = 0.036 and Ai2{{Z^A-'^/'^)/{A)) ~ 
-1.0. 

There are different possibilities for an evolving value of (j) to produce a 
variation in /i. If the variation is entirely due to the isospin-conserving coupling 
Xi then the differential acceleration vanishes, but in general one expects all 
couplings to be nonzero. (Variation of a also implies that A^ and Agm are 
nonvanishing.) We require (Ap — Ae)A(^ to be of order 10~^, where the change 
A(j) occurs between redshifts 2 to 3 and the present. If we have a nonlinear 
dependence of particle masses and a on (p then the functions Xe{4>) etc. must 
be integrated back to the value at redshift 2 to 3 and we obtain instead 
averaged values 

1 



Ae 



I Ap 



(where we earlier fixed 
(f) couplings we have 



0). In terms of the isospin conserving and violating 



Aln^ 



A,: 



1 



A^ - Ae A<p 



(8) 



If (p is not evolving fast, Acf) < 1 then at least one of these coupling constants 
should be of magnitude 10~^ or greater. 

If as in |38j the variation of fi is entirely due to the electron coupling then 
r] will be at least 10~^^. However, if both Aj and Ae are of order 10~^ then 
the lower bound on rj becomes much larger, a few times 10~^^. If Aj happens 
to vanish, then the value of rj will depend on the relation between A^,, Ae and 
Aem- The contribution of electromagnetic self-energy to the proton mass (which 
enters via A^,) is tiny and cannot lead to observable changes in fi, therefore Aem is 
bounded only by measurements of a. If we take a variation Aa/a = 0.6 x 10~^ 
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then the contribution of \em to r] is about an order of magnitude larger than 
that of Ae, comparing 0.036Ae/1837 with 7 x lO^'^Aem, where Ag is chosen to 
give A/i/yU of a few times 10^^. 



3 Scalar couplings in unified models 

In the previous discussion we ignored the underlying physics which should lead 
to relations between the scalar coupling strengths Aj, A^,, Ae and \em- In partic- 
ular the isospin- violating coupling A^ can only arise from the coupling of (p to 
up and down quark masses or electromagnetic self-energy. Now recall that the 
ratio niq/Mp, where Mp is the Planck mass, is mainly due to the hierarchically 
small value of the Higgs v.e.v. v/Mp, which is equally true of the electron mass: 
hence even without invoking unification it seems very likely that the functional 
dependence of rriu^d/Mp on (p will mirror that of nie/Mp. Symbolically 

Now the contribution of — to dninp is a few MeV, hence the value of A^ 
is suppressed by the ratio (m^ — rnu)/mN, relative to the quark couplings, and 
its contribution to r] will likely be comparable to that of Ag. 

To be more systematic, we note the following relations valid for small vari- 
ations in particle masses: 

Alng^ = Aln^ = £^A.A^. Aln^ = A.A^ (10) 

which will allow us to connect these low-energy quantities to underlying cou- 
plings of the cosmon. We choose to look at models with unification of gauge 
couplings in which the cosmon may couple to four quantities: the unification 
mass Mx, the unified fine structure constant ax defined at the scale Mx, 
the Higgs v.e.v. v, and (for supersymmetric theories) the soft super symmetry- 
breaking masses m, which enter as thresholds in the RG evolution equations. 
Their variations are defined as 

Mx 

A In — = AmA0, a In ax = Ax A0, 

t) — 771 — 

A\n— = XHAcP, Aln— = \sA<p (11) 
Mx Mx 



and in nonsupersymmetric theories we set A^ = 0. The results of |2E1 can now 
be used to relate the nucleon and electron masses to the underlying quantities. 
The electron mass varies as 

TTl — 

Alnj^ = Xm + Xh^(P (12) 

thus simply Ae = Xm + ^h] the same is true for the quark masses (neglecting the 
variation of Yukawa coupling constants). As an intermediate result we require 
the variation of the QCD scale Ac which satisfies 

Ah,^.A.A^=(|^-A, + ?A„ + lA.)A^ (13) 
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and the variation of a given by 



A In a EE \eml^4> - \ Ax + — 77 Aa + Ah + — As 

[ 6ax zvr \ 3 3 

Then the averaged nucleon mass coupling is 

\ = Aa/ + ^Ax + (0.39 ± 0.12)Ai^ + 0.35As 

ax 



A0. 



(14) 



(15) 



where hght quark mass contributions, heavy quark and superpartner thresholds 
are accounted for; the uncertainty in the Xu term arises from the strangeness 
content (see for example 00]). The isospin violating coupling is found from the 
expression 



np 



2.05 1 + Aln 



Mp 



0.76 ( 1 + Alna + Aln^ 



with coefficients derived from this leads to 

A„ ~ 1.4 X 10-3 [\m - 10.2Ax + IA&Xh - O.27A5] 



MeV 
(16) 

(17) 



where we take ~ 24. Note that that the contribution of varying a to 
this mass splitting is always very small compared to that of Ac which is much 
more sensitive to the underlying parameters. We see immediately that the 
contribution of A„ to the variation of /i is negligible; we have 



A In /i ~ [13.0Ax - (0.61 ± {).12)\h + 0.35Acj] A(/.. 



(18) 



We can now for a range of unified models determine the relations between 
variations of ^ and a and the low-energy coupling constants that give rise to 
differential accelerations. We ignore \m since this coupling does not affect 
either a or /i. 

Scenario 1: varying unified coupling We first consider the case when only 
Ax is nonvanishing, among the underlying cosmon couplings. Then we have 



A In a ~ 0.52Ax A(/), A In // ~ 13Ax A(/), 



(19) 



Ae = 0, Ai = 13Ax, A^ = (-1.4 X 10-2)Ax (20) 

thus A In /I ~ 25A In a (which would put the variation of a beyond the reach 
of current observations). The magnitude of Ax is at least of order 10~®, to 
produce a nonzero variation of /i. Now the differential acceleration arises from 
the product XiXy/2 which is at least 10"^^, leading to rj of magnitude just below 
10"^^; the contribution proportional to Xem is slightly smaller. 
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Scenario 2: varying Higgs v.e.v. without SUSY Next we consider if the 
variation were to arises solely from a change in the Higgs v.e.v. relative to the 
unified scale, thus Xh is nonzero. Then we have 

Alna ~ (2 X 10~^)AhA^, A In^i ~ -O.GA^Ac^, (21) 

Ae = Xh, Xi ~ OAXh, Xy = 2x IQ-^Xh (22) 

In this case the fractional variation of a should be about 300 times smaller 
than that of /x! Now to reproduce a nonzero variation of /i we require Xh to 
be a few times 10~^ and AjA^/2 is again at the level of 10~^^, with the electron 
contribution AjAg/ (2 • 1837) only slightly smaller. Again r] is bounded below by 
lO"^"^, given a nonzero variation of fi. 

Scenario 3: variation of SUSY-breaking hidden sector and Higgs 

v.e.v. Next we consider what happens if supersymmetry-breaking arises from 
a hidden sector with a mass scale that is generated by gauge dynamics, and 
assume that electroweak symmetry breaking is triggered by soft SUSY-breaking 
masses. We will first take the hidden sector gauge coupling to vary indepen- 
dently of the visible sector gauge couplings, thus 

Ax = 0, Xs^Xh. (23) 

Then we find 

Aln^ = (-0.3±0.1)AsA(^, Alna = (1.3 x lO"2)AsA0, (24) 

Ae~A5, A, ~0.74As, A, = 1.7 x lO-^A^, (25) 

with As being at least 10~^ to produce a variation of fi at the observed level. 
Here the fractional variation of a should be about 20 times smaller than that of 
/J, and of opposite sign. Then the product AjAt,/2 is at least a few times 10~^^ 
and r] must reach the level of 10~^^ level, pushing current observational limits. 
The Ae and Xem contributions are slightly smaller. 

Scenario 4: varying unified coupUng plus SUSY-breaking Finally we 
consider what may happen if the hidden sector gauge coupling varies with the 
unified coupling at the fundamental scale. As first pointed out in |25j this leads 
to the relations 

As-A//~34Ax. (26) 

This results |28j in Aln/u being of the same magnitude as A In a. More con- 
cretely 

Aln^ = (4±4)AxA(^, Alna~l.OAxA0 (27) 

thus Ax is just below 10^^, or larger, to match the claimed nonzero variations. 
We then find 

Ae~34Ax, Ai~38Ax, A„ ~ (4 x 10-2)Ax, (28) 
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thus the products of couphngs appearing in r] are a few times 10~ , leading 
to differential accelerations of order 10~^^, already outside the experimental 
range! 

This comes about because the variations of nip/Mp and mg/Mp are large 
but cancel against one another in /i; however the differential accelerations due 
to the nucleon mass splitting, the electron mass and the nuclear electromagnetic 
energy do not in general cancel against each other. The contributions from Ae 
and Aem are slightly smaller than those from A^, and may have opposite sign, 
but even if they did cancel in one Eotvos-type experiment this could not hold in 
other setups where source and test bodies had different values of Z and A. This 
scenario could only survive if the variations of a and /i were below currently 
claimed nonzero values. 

One can also consider a non-SUSY GUT scenario where the unification mass 
scale and coupling are Mx — 10^^ GeV and ~ 40, where the electroweak 
scale would arise from some technicolor- like sector, giving A^ = (31 it 2) Ax if 
the technicolor group is unified with the Standard Model. The results for this 
type of model are in fact very close to Scenario 3. There is a cancellation in 
the (fractional) variation of /i such that is comparable to the variation of a, 
however the scalar couplings are relatively large and rj is of order 10~^^. 

Models where there is a partial cancellation between variations of mg/Mx 
and rUp/Mx have been put forward in and the last reference of as 
a way of reconciling the claimed nonzero variations of a and ^ which are in- 
compatible in many simple unified theories. However, these estimates of WEP 
violation show that such models are likely to have problems to respect current 
experimental bounds. 



Simplified expression for WEP vs. varying mu We have seen that in all 
these scenarios based on unification the source factor (first bracket) in Eq. ((Tj) 
is strongly dominated by the isospin-conserving scalar coupling Aj, while the 
second factor depends on the electron coupling, the isospin-violating coupling 
A^ and the electromagnetic coupling Xem- The isospin-violating contribution to 
A In /i in Eq. Q is negligible. Therefore the relation between the variation of 
jji and the differential acceleration rj can be simplified to 



J? > 



(Aln/i)2 ~ 2{Xi-X,f 



A 



1837 



/72 /i-l/3\ 

Ai2/p + 7x 10-%„Ai2-^ ^ 



{A) J 
(29) 

due to the fact that all scalar couplings arise from one underlying parameter 
whose variation, for any given Aln^u, is proportional to (A, — Ag)"^. 



4 Evolution of the cosmon or quintessence field 

Since the estimates of the previous section are interestingly close to the current 
experimental limit of r], and we made only rough estimates using the assumption 
that A(j) was less than unity, we will refine them by examining in more detail 
the cosmological evolution of (p and applying constraints on its kinetic energy 
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density. In previous work j21j a physical model was adopted with an exponential 
potential and a general (non-canonical) kinetic term allowing the past evolution 
of the field to be found in terms of a few parameters. However in general the 
evolution and variation of "constants" may depend on undetermined functions 
of a scalar field, either in the scalar action (kinetic terms plus potential) or its 
couplings to matter. 

Another question is whether the evolution is homogeneous, that is, whether 
the scalar value depends only on cosmological time, or is position- or environment- 
dependent. This question was investigated for gravitationally collapsing regions 
of spacetime in |32| and more recently for virialized structures in , where the 
result was found that the rate of change of (p is the same within virialized regions 
as in the background and the value of the field depends only weakly on position, 
a result also obtained in |341l35j (confirming the analysis of |^). However, note 
that if a scalar field couples strongly to matter such that its expectation value 
and mass are determined by the local density, both the cosmological evolution 
and the possible scalar-mediated forces behave quite differently |36| \'A7\ . It has 
been claimed that such a "chameleon" field cannot give rise to a fractional vari- 
ation of a at the 10~^ level, but that bound may be due to particular choices 
made in setting up the model of ^331, rather than a general result. We will 
consider here only the weak coupling regime where the mass of the scalar field 
is determined primarily by its potential V{(j)). 

4.1 Relating WEP violation to atomic clock measurements 

The (j) couplings to matter can be related most directly to the present rate 
of change of dimensionless couplings, which we now derive. For a canonically 
normalized scalar field we have kinetic energy T = {\/2)(j)^ and a potential 
V{4>), resulting in the energy density = T + V and equation of state = 
{T — V)/{T + V). There may also be contributions to the scalar equation 
of motion from interactions with matter, dark matter etc. and as shown for 
example in [36^ these may alter the effective equation of state to be less negative. 
However we do not expect such contributions to greatly affect our bound, which 
is derived from the maximum kinetic energy of </> consistent with observation 
and does not depend strongly on w^. 

The energy density fraction Q.^ is defined as p^j Pc = P(^/3-?^^Mp. Hence we 
obtain 

|: = ^317^(1+^. (30) 

The usefulness of this equation is that the Hubble constant is now well measured 
whilst the quantities inside the square root on the RHS are bounded above. 
Given current cosmological observations, there is an allowed region in the O^- 
plane inside which Eq. (|3U() attains a maximum value, which we can estimate 
roughly as ^3 x 0.75 x 0.2 ~ 0.7. Using ~ 7 x 10"^^y"^ we find 

^ < <^max =^ 5 X 10-"y-\ (31) 

For a given unified model where all varying quantities are determined by a 
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single parameter A we can write 



(y. — jj, — 

- < ClA(/)max, - < C2A(/)max, (32) 

a fi 

where the numerical constants ci and C2 arise from relations detailed in the 
previous section. Similarly the results for rj can be summarized in the form 

^ = K^^X" (33) 

for some constant K, where for simplicity we omit the contribution of nuclear 
electromagnetic energy which was subdominant in all the scenarios considered. ^ 
Thus we find 

^ - 20Lx '2 W) ' 4 V3.7 X ^'^^ 

taking Ai2/p to be 0.036 as before. An exactly analogous inequality follows by 
replacing with a and C2 with ci. 

Now if we suppose that a or /i are nonzero at present (either due to direct 
measurement or to a particular model of cosmological evolution) then in any 
given scenario of unification we find the values of K/cl and i^/c^ and obtain 
the lower bound on r]. 

Conversely, given bounds on rj, for any given unification scenario we can 
obtain bounds on a/ a and fi/fi which must be satisfied for a cosmological 
variation arising from a scalar field effective action. These may be directly 
compared with atomic clock measurements. Note that a nonzero result for rj, 
while extremely interesting, would not have direct implications for "varying 
constants" within this framework since the rate of change of (p is not bounded 
below. 

As an example we evaluate such bounds within Scenario 1, where the rel- 
evant scalar coupling is Ax and the numerical constants are K/c1 ~ 0.67 and 
K/c2 — 0.0011. Thus for a hypothetical variation a/a ~ 10~^^y~^, we find 
that rj should be at least 5 x lO"^"*. 

For Scenario 4, which is barely ruled out from the point of view of producing 
the claimed nonzero variations of a and fi, we find the larger values K/cf ~ 1.5 
and K/c2 — 0.1 reflecting the fact that such scenarios produce small /u variations 
through partial cancellation. 

4.2 Bound on variation of (f) at earlier times 

Equation (|3n)l may also be extended to earlier times, given bounds on the energy 
content of the universe derived from cosmological measurements. Changing 
variables from time to redshift we find for the change in the scalar field value 
back to some redshift z 

|A0| < ln(l + z) max (^^Sfl^l + w^)^ (35) 



^With an experiment using test bodies of very large Z, this contribution could become 
competitive or dominant. 



11 



where the maximum is to be evaluated over the range from z to the present. 
A very conservative bound is a scalar field behaving like cold dark matter, for 
which = and 17 < 0.35, thus the square root is 1 or less. This leads to a 
limit such as |Ai^| < 0.12 at the Oklo epoch z = 0.13 (see the last reference 
of [S]), or |Ai^| < 1.4 at redshift 3. This justifies the previous assumption that 
|A0| was at most of order 1 at redshifts relevant for QSO absorption spectra. 
We emphasize that in almost all viable scalar theories such values of Ac/) will 
be a considerable overestimate, since it would require considerable fine tuning 
to saturate these upper bounds on <^ at all times. 

In theories where the nonlinear dependence of coupling strengths on the 
canonically normalized scalar (\) is known, such results can be used to check 
consistency of measurements of a, [i and other quantities at different epochs. 
In general as stated earlier the dependence is allowed to be an arbitrary function 
of and only the derivative at the current epoch is probed by WEP and atomic 
clocks. 

5 Conclusions 

We have investigated bounds arising from violations of the Weak Equivalence 
Principle in scalar field theories in which the proton-electron mass ratio /i may 
vary. We related such bounds to unified theories, astrophysical observations of 
quasar absorption spectra and atomic clock measurements of the constancy of 
physical parameters. Model-independent bounds on the rate of change of the 
scalar field and on its total variation at any given epoch are obtained from the 
effect of its kinetic energy on the expansion of the Universe. 

For a given nonzero variation of /x, the scalar couplings to the proton or 
electron mass are bounded below, thus so is the resulting differential accelera- 
tion, up to a single parameter Kjc^. This parameter is a ratio of scalar coupling 
constants, which arises from how the nucleon and electron masses are related 
in the underlying particle theory. Hence such underlying theories, in particular 
gauge unification at high energy, could be tested with variation of ii. 

Since the value of a may also vary in such underlying theories, we also ana- 
lyzed its variation in relation to WEP violation, extending previous works. We 
found that in many scenarios for generating varying /i from unified models, the 
differential acceleration parameter r/ is near or above its current experimental 
limit. In particular, in scenarios where due to cancellations the variation of \i 
is not much larger than that of a, the contributions to differential accelerations 
do not cancel and the violation of WEP is at or above the experimental limit. 
Thus it seems very difficult to reconcile currently claimed nonzero variations of 
/X and a within unified models. 

Our results differ from those of |HSI, where the cosmological variation of \i 
was bounded to be many orders of magnitude below current levels of sensitivity. 
The main difference is that in I55v the source of the scalar field variation was 
taken to be the coupling to the local electron density, whereas we allow the 
field to be driven by a potential energy and in principle by couplings to other 
constituents of the Universe. This illustrates the point that in models where the 
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form of the scalar Lagrangian is known or assumed, more stringent bounds may 
be found. Our aim was however to find bounds which are as model- independent 
as possible. 
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